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^ ■ Abstract 

We review the pion-photon transition form factor calculated in the Schwinger- 
I Dyson approach and an impulse approximation. We present results for it far 

ly^ I above the scales presently accessible for measurement, up to 36 GeV^, and 

' demonstrate agreement with the analytically inferred asymptotic behavior, 

Q>^ ■ for which we also provide a new derivation. We discuss how measurements at 

Jefferson Lab can provide information on how quarks are dynamically dressed. 

1. INTRODUCTION AND SURVEY 

^ I A modern version of the constituent quark model with many attractive features is one of 

• remarkable achievements of the Schwinger-Dyson (SD) approach (reviewed in Refs. [0,0) to 
the physics of quarks and hadrons. It is well-known that most quark-antiquark bound state 
approaches have grave problems when they are faced with the electromagnetic processes 
dominated by Abelian axial anomaly. (See Ref. for a brief review and more references.) 
However, these problems are resolved in the SD approach thanks to its good chiral properties. 
In the chiral limit, the axial-anomaly result for 7r°(p) —>■ 'y{k)'-f{k') transition amplitude 

T,.(0.0) = j^ (1) 

is reproduced analytically and exactly in this approach [|,||. The amplitude in Eq. (D) is 
the limit when both photons are real (A;^ = k''^ = 0), of the general amplitude TT^o{k'^, k'"^). 

At first, 7r° ^^(*) processes and other applications were studied in the approach with 
physically motivated Ansdtze {e.g. 0,|]), instead of the solutions of SD equations (SDE) with 
specified dynamics. However, the exact evaluation of Abelian axial anomaly in the closed 
form (|1]) occurs also in the variant where one actually solves SDE for the quark propagators 
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and, in a consistent approximation and with the same model interaction, Bethe-Salpeter 
(BS) equations (BSE) for quark-antiquark bound states. This is the consistently coupled 
SD-BS approach, of which [0 and references therein provide an elaborated example. 

Since the anomaly is independent of hadronic structure, it is essential to note that the 
relation (|1|) is successfully reproduced in the SD approach irrespective of what concrete 
Ansatz for the quark propagator 

S{q) = [A{q'U-B{q')]-\ (2) 



is used, or what concrete consistently coupled SD-BS solutions are employed for the dressed 
quark propagator (|) and the pion qq bound-state BS vertex Tj^o. 

The successful treatment of the Abelian axial anomaly is possible in the SD approach 
because this approach incorporates the dynamical chiral symmetry breaking (D^SB) into 
the bound states consistently, so that the pion, although constructed as a qq bound state, 
appears as a Goldstone boson in the chiral limit. With that, the bound-state descriptions 
of mesons are reconciled with chiral requirements stemming both from QCD as underlying 
fundamental theory, and from phenomenology. The scenario is essentially of the Nambu- 
Jona-Lasinio (NJL) type, but without its low cutoff. Consequently, the quark-antiquark 
bound state mesons can finally provide adequate descriptions of the processes where the 
axial anomaly is important {e.g., see [^H^,3~|l0|,|3|| ) , such as the two-photon processes of 
light pseudoscalars, of which the tt^ —>■ '-ff decay is the cleanest example of an anomalous 
electromagnetic decay. 

Having gotten in hands the fully satisfactory result (|1]) for the 7r° —>■ 77, the next thing 
to explore was its one-off-shell-photon extension T^o(/c^, 0) in the SD approach. Frank et al. 
@] studied it using Ansdtze for dressed quark propagators (0), and then Ref. did it in 
the consistently coupled SD-BS approach, utilizing the solutions of Ref. [0]. However, both 
papers called for additional studies in two respects. First, both had employed the soft and 
chiral limit for the pion; i.e., the BS vertex was approximated by its leading, piece: 

r.o(g,p)^r.o(g,0)=75A3^^ • (3) 

Jtt 



(A^ is the third Gell-Mann matrix of flavor SU{3).) In addition, Refs. could compare 
their results just with then only available data by CELLO at < 2.5 GeV^. 

Lately, however, the interest in the form factor Tn-o(— (5^,0) for the transition 
'y*{k)'~f{k') — i> n^{p) (where fc^ = —Q^ 7^ is the momentum- squared of the spacelike off-shell 
photon 7*), has again been growing for both experimental (the new CLEO data [|1^] and 
plans for new TJNAF measurements and theoretical reasons - e.g., see ||14|| . 

At large values of Q^, the 7*7 7r° transition form factor should be adequately given 



by perturbative QCD (pQCD) - e.g. see [p!5HT7[]. Nevertheless, it is still not quite certain 



which Q is sufficiently large. For example, according to Refs. []TB|, just pQCD may still 
be not quite sufficient even at the highest of the presently accessible momenta, < 10 
GeV^. The pQCD approaches start having problems as decreases more and more into 
the nonperturbative domain. Also, see Sec. IV of the recent Ref. fl^ for clarifications how 



such approaches |]T6| fail to reproduce the anomaly-induced value (|I]) at = 0. 
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On the other hand, the treatment of the axial anomaly is a strong point of the SD 
approach, as pointed out above. What must then be explored in the coupled SD-BS approach 
is if the large behavior is satisfactory. In particular, the comparison with the new 
CLEO data - at up to 8 GeV^ - must be made. Whether the large behavior 
is satisfactorily close to the data and to the predictions of pQCD is a tricky question for 
a constituent quark model. Namely, the calculation of the transition form factor carried 
out in the simple constituent quark model (with the constant light-quark mass parameter 
m^), leads to Tjro(— (5^,0) oc (m^/Q^) ln^(Q^/m^) as — >• cxo, which overshoots both the 
CLEO data and pQCD predictions considerably This is because of the additional 

ln^((5^)-dependence on top of the large behavior 

7;o(-Q^ 0)=J^ iJ^ const as ^ oo), (4) 



favored by pQCD and other QCD-based theoretical predictions such as [^5]-|20| and references 



therein. Publication of the CLEO data |T2| made it clear that the large-Q behavior is 
also favored experimentally |12 . 

Fortunately, the constituent quark model provided by the SD approach, in which 
B[q^) / A[q^) plays the role of the dynamically generated g^-dependent mass, has turned 
out not to suffer from that shortcoming. Namely, we have shown in Ref. [^] that the above 
asymptotic behavior (|) is in this approach obtained in the model-independent way. 

In this paper, we provide yet another way to obtain the asymptotic — > cxd behavior 
found in Ref. Predictions for T.^o^—Q'^^Q) for finite values of can be obtained 



by picking a definite model for interactions between quarks and obtaining corresponding 
solutions of the SDE and BSE. In this paper, just as we did in Ref. [^, we use the model 
of Ref. 0. However, here we give the results over much wider range of - see Fig. 1 
and the discussion thereof. Comparison with the data reveals in what way more precise 
measurements at intermediate momenta, feasible at TJNAF, can both give insights in the 
hadronic structure and provide guidance how to improve presently existing models in the 
SD-BS approach. In addition, the coupled SD-BS variant of the SD approach (such as that 
in Ref. [^) can readily be extended and applied beyond the soft and chiral limits [given by 
Eq. (^] for the bound-state vertex - e.g. see Refs. P, |2^ , p!0| j3[] . We did this when calculating 
^7ro(— <5^, 0) in Ref. |2l|], but here we locate and point out the main reason for the difference 
between the full calculation of the transition form factor and the one using the the soft and 
chiral limit approximation (^ for the pion BS vertex. 



2. DRESSED QUARKS AND BS- VERTICES NEED DRESSED qq-f VERTICES 

In the coupled SD-BS approach, the BSE for the pion bound-state qq vertex r^o(g,p) 
employs the dressed quark propagator (0) obtained by solving its SDE. Solving the SDE and 
BSE in a consistent approximation is crucial {e.g., see Refs. P3| , pi| ,[7|, P3| , P^ ) for obtaining qq 
bound states which are, in the case of light pseudoscalar mesons, simultaneously also the 
(pseudo-)Goldstone bosons of D^SB. 

Following Jain and Munczek p3|j2^ j7[|, we adopt the ladder- type approximation some- 
times called the improved ||2^ or generalized ladder approximation (employing bare quark- 
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gluon-quark vertices but dressed propagators). For the gluon propagator we use an effec- 
tive, {partially) modeled one in Landau-gauge [^^ [7|, given by G{—P){g'^^ — IH^ /t^) . 
(This Ansatz is often called the "Abelian approximation" [^.) The effective propaga- 
tor function G is the sum of the perturbative contribution Guv and the nonperturbative 
contribution Gir. G{Q'^) = GuviQ"^) + Gir^Q"^) , (Q^ = -/^) • The perturbative part 
Guv = (167r/3)as(Q^)/Q^ is well-known from perturbative QCD, so it is not modeled 
2|,|2|J3jg|Tgg]. As in Refs. B|TU|,§, we follow Refs. [^^0 and employ the two-loop 



asymptotic expression for ^^(Q^). For the modeled, IR part of the gluon propagator, we 
adopt from Ref. [0] Gir{Q'^) = (167rV3) a Q^g-z^Q', with their parameters a = (0.387 GeV)"^ 
and = (0.510 GeVy^. For details of how we solve SDK and BSE, we refer to Refs. [|,|T0|,|. 
To high accuracy, we reproduce Jain and Munczek's solutions of SDE for the dressed 
propagators S{q), i.e., the functions A(g^) and B{q^), as well as the BS solutions for the 
four functions comprising the pion bound-state vertex F^o. 

The 7r°7*7 transition tensor {T'^o) and scalar (T^o) amplitudes, related by 

T^:^{k, k') = e^^^^k^k'pTAk\ k") , (5) 

are found in the present paper as in Refs. PJTO|J^. That is, the pseudoscalar- vector- vector 
(PVV) triangle graph is calculated by employing the framework advocated by (for example) 
HH6|,^,|8| in the context of electromagnetic interactions of BS bound states, and often called 
the generahzed impulse approximation (GIA) - e.g., by P,|27[]. To evaluate the triangle graph, 
we therefore use the dressed quaxk propagator S{q) and the pseudoscalar (P = 7r°, rjs, rjo, rjc...) 
meson BS bound-state vertex Fp(g,p). (See Fig. 1 in Ref. |]2T[], for example.) Equivalently, 
we can work with the BS-amplitude Xpil^p) = S{q + p/2)V p{q,p)S{q — p/2). This is in fact 
what we do, following Ref. 0. 

The third ingredient crucial for GIA's ability to reproduce the correct Abelian anomaly 
result (^, is employing in the triangle graph appropriately dressed electromagnetic vec- 
tor vertices r^{q',q), which satisfy the vector Ward-Takahashi identity (WTI) {q' — 
q)^T^{q',q) = S~^{q') — S~^{q) . Namely, assuming that photons couple to quarks through 
the bare vertex would be inconsistent with our quark propagator S{q), which, dynamically 
dressed through its SD-equation, contains the momentum-dependent functions A(g^) and 
B[q'^). The bare vertex 7^^ obviously violates the vector WTI, implying the non-conservation 
of the electromagnetic vector current and charge. Solving the pertinent SD equation for the 
dressed quark-photon-quark [qqj) vertex F'^ is still an unsolved problem, and using the 
realistic Ansatze for F^ is presently the only practical way to satisfy the WTI. Following 
Refs. [i|7|,|,i (for example), we can choose the Ball— Chiu (BC) [|28|| vertex; i.e., — ^bci 



= A^q'^i')^ + ^^^,{A.{q'\q')^l±ll- - B_{q'\q')} , (6) 



where H±{q''^ ,q^) = [if(g'^) ± if(g^)], for if = A or B. This particular solution of the 
vector WTI reduces to the bare vertex in the free-field limit as must be in perturbation 
theory, has the same transformation properties under Lorentz transformations and charge 
conjugation as the bare vertex, and has no kinematic singularities. It does not introduce 
any new parameters as it is completely determined by the dressed quark propagator S{q). 

Another WTI-preserving choice for F'^ can be a vertex of the Curtis-Pennington (CP) 
type, i.e., F^ = V^^. + AF'^ where the the transverse addition AF'^ is of the type 
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2d{q',q) 

Two especially suitable Ansdtze for the dynamical function d{q', q), ensuring multiplicative 
renormalizability of fermion SDE beyond the ladder approximation in QED4, are 



d^{q',q) = ^^J{q"±qr + 



B\q'^) B\q^ 



(8) 



The original CP Ansatz = F^p employed d{q',q) = d-{q',q) p9|. We will use it in 
analytic calculations of T{—Q^,0), which are possible for = and — > 00. However, 
in the numerical calculations, which are necessary for finite values of 7^ 0, we prefer 
to restrict ourselves (besides the minimal BC vertex) to the modified CP (mCP) vertex, 
Ti^QP, resulting from the choice d = d+ in Eq. (0). Namely, as pointed out in Ref. the 
numerical calculation of T(— Q^,0) employing F^^^-p vertices is free from certain numerical 
difficulties arising, in this application, from the CP denominator function d-{q',q). 

In contrast to the BC one, the mCP vertex is also consistent with multiplicative renor- 
malizability, like the original CP vertex. In the present context, the important qualitative 
difference between the BC-vertex on one side, and the CP vertex as well as the modified, 
mCP vertex on the other side, will be that T^(j{q' , q) 7^^ when both g'^, q^ ±00, whereas 
^cpW 1 1) ~^ 1^ F(^^p(g', q) 7^ as soon as one of the squared momenta tends to in- 
finity. This turned out to lead to the same coefficient of the Q"^ 00 behavior for 
the CP vertices and mCP vertices, namely = 4/3, but a larger one for the BC-vertices. 

We have checked that the = case, i.e., the vr" — > 77 amplitude (|I]), is reproduced 
analytically employing the CP vertices and mCP vertices in the same way as when the BC 
vertices were employed in earlier applications, e.g. P,|6|,|9|JT0| . 

In the case of 7r°, CIA yields {e.g., see Eq. (24) in Ref. [|T0|) the amplitude T^oik, k'): 



T^,-{k, k') = -N,^ f -^tiiTf'iq -^,k + q- ^)S{k + q-^ 
' ' 3V2J (27r)4 ^ 2' ^ 2' ^ ^2 

p p 



X F^(fc + g - |, g + |)x(g,p)} + {k^k',fl^u). (9) 

Here, x is the BS amplitude of both uu and dd pseudoscalar bound states: x = Xuu = XdS 
thanks to the isospin symmetry assumed here. This symmetry likewise enables us to continue 
suppressing fiavor labels also on the quark propagators 5" and qq^y vertices F^. We follow 
the conventions of Ref. |jlO[, including those for the fiavor factors and fiavor matrices A". 
Then, XTTo{q,p) = x{q^p) / V2, so that the prefactor l/3\/2 in Eq. (H) is just the fiavor 
trace tT^Q'^X^/V^) where Q = \^/2 + \^/2^/Z = diag{+2/3, -1/3, -1/3) is the quark charge 
matrix. 



3. OUR RESULTS AND COMPARISON WITH OTHERS 

Already in Ref. [^] (where only the BC vertex was employed), the transition form factor 
T^o(— (5^,0) was numerically evaluated (for < < 2.8 GeV^ only) employing - for the 
pion only - the soft and chiral limit (0). 
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In the subsequent work ||2T|, we went beyond this approximation, using our complete 
solution for the BS vertex r^o(g,p), viz. the BS amplitude X-rTo{q,p), given by the decom- 
position into 4 scalar functions multiplying independent spinor structures. We calculated 
the transition form factor not only for the BC vertex, but also for the mCP vertex. The 
approximation we kept in Ref. |21| (and which we still keep in the present paper), is dis- 
carding the second and higher derivatives in the momentum expansions of the SD solutions 
Ai^q"^) and -B(g^) and BS solutions xil^p)- R-sf. Tt,q{—Q'^,0) was evaluated up to 
Q"^ ^ 8 GeV^, which is roughly the limit of the presently accessible transferred momenta. 
In this paper, we evaluate the transition form factor up to 36 GeV^. We present curves for 
(32T^o(-Q^0) evaluated in various ways first in Fig. 1, depicting their evolution all the way 
to 36 GeV^. This squeezes the presently existing experimental points into the first quarter 
of the scale; however, the issue of comparing and improving the agreement with the data at 
the accessible momenta will be returned to later, and now we want first to illustrate how 
various theoretical transition form factors at finite connect to the asymptotic behavior. 

In Fig. 1, the sohd line is our prediction for Q^T^o{—Q^,0) evaluated with the BC 
vertex. It reaches the asymptotic behavior faster than other curves. Already after ^ 4 
GeV^, the tiny changes cannot be observed in the solid line, but only in the corresponding 
numerical data, as the asymptotic value is approached from below. For comparison, the 
dash-dotted line gives the results obtained (also with the BC vertex) in the approximation 
(§) of soft and chiral limit as explained above. At ~ 4 GeV^, where the Q^T^o(— Q^, 0) 
from the full calculation already practically reaches the asymptotic limit, the difference 
between it and Q'^T^o^^{—Q'^,0) from the approximation (|^), is somewhat above 10%. This 
difference gets bigger as grows, and asymptotically it reaches some 22% of the limit value 
of Q^T^o{—Q'^,0). However, this still means that the effect of the soft limit approximation 
on the transition form factor is much less than on the charge form factor. Namely, Maris and 
Roberts pointed out that omitting the pseudovector components of the pion lead even 
to the wrong asymptotic behavior, as instead of l/Q"^, in the case of the charge form 
factor of the pion. In the present case of the transition form factor, the correct l/Q^ leading 
behavior is nevertheless obtained not only for the results of the full calculation, but also in 
the soft limit approximation (|^). It is just that the coefficient of 1/Q^ is underestimated 
with respect to the full calculation. 

The dashed line was obtained in the same way as the solid line, but employing the mCP 
vertex. Even at = 36 GeV^, the asymptotic behavior is obviously not yet reached in 
the case of the mCP vertex. However, later we will comment on how Ref. 11211 showed 



analytically that the behavior given by Eq. (^) must be reached at some point, although at 
much higher values of when the mCP vertex is used. 

How do our transition form factors agree with other theoretical approaches? The vector- 
meson dominance (VMD) model and the QCD sum rule approach [|l8l give the transition 
form factor which is some 10% below our "BC" T^o(— Q^, 0) for the presently largest accessi- 
ble values of Q^, i.e., around 8 GeV^. Therefore, in that region, our BC- results are between 
the uppermost line in Fig. 1 (the dotted constant line at 2/^) denoting the asymptotic pQCD 
l&i version {J = 2) of Eq. (I), and the results of VMD (e.^., see fll), the recent pQCD 



calculation by Ref. as well as the QCD sum rule results of Radyushkin and Ruskov [jl8| 
Q'^Tt^o {—Q"^ , 0) of VMD and pQCD rise with Q^, albeit with different rates, while that from 
the QCD sum rules of Ref. |T^ starts almost imperceptibly falling after ~ 7 GeV^. 
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The large-Q leading power-law behavior (m was first derived from the parton pic- 



ture in the infinite momentum frame - e.g., see Ref. In this and other similar 

pQCD approaches, the precise value of the coefficient of the leading l/Q"^ term depends 
on the pion distribution amplitude (pTr{x) which should contain the necessary nonpertur- 
bative information about the probability that a partonic quark carries the fraction x of 
the total longitudinal momentum. E.g., the well-known example of a "broad" distribution 



= f-,r5V3x{l — x){l — 2xY (proposed by Chernyak and Zhitnitsky motivated by 
sum-rule considerations) leads to J' = 10/3, but it is too large in the light of the latest data 
T^. In contrast, the asymptotic distribution (p^{x) = f-^^/Sx^l — x) (depicted by dotted line 



in Fig. 2 and favored by Lepage and Brodsky |T^) yields J = 2, resulting in the constant 
dotted line crossing the upper error bar of the presently highest-Q^ data point in Fig. 1. 

In the strict\\y{Cp'^ oo limit, every distribution amplitude must evolve into the asymp- 
totic one, '^■„[x) Lp'^{x), if the effects of the pQCD evolution are taken into account [jl5 . 
However, even at Q^-values larger than the presently accessible ones, other effects may still 
be more important than the effects of the pQCD evolution. This is the reason why other ap- 
proaches and other forms of '^t^[x) should be considered even when they do not incorporate 
the pQCD evolution. This is in line with Radyushkin and Ruskov's ( |jl8| and their refer- 
ences therein) pointing out desirability of having direct calculations of Tt^o{—Q'^, 0) without 
a priori assumptions about the pion distribution amplitude y^^. One can then consider the 
opposite procedure from the one which is standard in pQCD: from such direct calculations 
of T^o(— (5^,0), one can draw conclusions about the distribution amplitude ip-,,. 

The form 

is suitable for representing various distribution amplitudes because it is relatively general 
rB|: C > 1 yields the distributions that are "peaked" or "narrowed" with respect to the 



asymptotic one (^ = 1), whereas C < 1 gives the "broadened" distributions which, however, 
now seems to be ruled out for the same reason as ip^^{x) quoted above, since JT" > 2 is ruled 
out empirically by CLEO ||12|. Namely, it is easy to see that Eq. (p!0[) implies ( = 2/ (3^7 — 4). 



In Fig. 2, we plot the distribution amplitudes of the form ([ToD for the three cases that 
are the most interesting for the present discussion. As seen below, these cases correspond 
to the C-values equal to 1 (dotted line), 1.5 (solid line), and 2.5 (dash-dotted line). 

At the highest presently accessible momenta, the asymptotic prediction {J' = 2) is 
lowered by some 20% by the lowest order QCD radiative corrections |31|, amounting to 



J' ^ 1.6, which fits the CLEO data well. Of course, these QCD corrections mean that 
(5^T^o(— Q^, 0) is not strictly constant, but according to Ref. ||3ll| it rises towards 2/^ so 
slowly that we can take it constant in practice. [Eqs. (19) and (18) from Ref. |3l|] lead to 
the line of empty squares in Fig. 1, that is, to (5^T^o(— Q^, 0) which grows just 4% from 
= 9 GeV^ to = 36 GeV^ when the number of relevant flavors is set to Nf = 5. 
The full calculation which employs the momentum- dependent Nf |^^, seems to give 
even slower variation of Q'^T^^o - see their Fig. 2.] Similarly slow variation results from the 
sum-rule approach of Radyushkin and Ruskov |T8|, which yields the transition form factor 
which is, on the interval from 3 to 8 GeV^, quite close to those of the pQCD approach of 
Brodsky et al. [^. The sum-rule |]TB[ Q^T^o(— Q^, 0) starts actually falling after ~ 7 
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GeV^, but so slowly that Eq. with the constant JT" ~ 1.6 represents it accurately at the 
presently accessible values of . This, just like roughly the same J associated to Ref. 



corresponds to a rather narrow distribution ( jlOD with C, ~ 2.5, depicted by the dash-dotted 
line in Fig. 3. 

The leading large-Q^ behavior as in Eq. (^), was obtained also by Manohar |^ using 
the operator product expansion (OPE). According to his OPE calculation, the coefficient 
in Eq. (^) giving the leading term is J = 4/3, which is below our large-Q^ T^o(— Q^, 0) by 
~ 20% when we use the BC vertex, but exactly coincides with the Q"^ ^ oo limit obtained 
when we use the CP or the mCP vertex - see Ref. |^T| and comments below. The coefficient 
J = 4/3 is the lowest one still consistent with the form (|T0|) because it corresponds to 
C = oo. The pion distribution amplitude (|10|) then becomes infinitely peaked delta function: 
M^) = (/./2V3)(5(x-l/2). 

For > 4 GeV^, our "BC" T^o(— (5^,0) also behaves in excellent approximation as @, 
with J' ^ 1.78. This would in the pQCD factorization approach correspond to (pn{x) (p^ 
with C ^ 1.5. On the other hand, our "mCP" T^o(— Q^, 0) falls off faster than for even 
the largest of the values depicted in Fig. 2. However, it does not fall off much faster, 
as our "mCP" Q^T^o{-Q^,0) at = 18 GeV^ is only 6%, and at the huge = 36 GeV^ 
is only 10% smaller than at = 9 GeV^ (roughly the highest presently accessible Q'^). 
Moreover, Ref. |21] showed analytically that, generally, the l/Q^-behavior of Eq. (^ is at 
some point reached in our approach, although Fig. 1 shows that for the mCP qq'y vertices 
this can happen only at significantly higher than it happens for the BC qq'y vertices. 



DISCUSSION, COMMENTS AND CONCLUSIONS 



The most important novelty in Ref. is its result on the large-Q^ asymptotic behavior. 
Namely, it implies that the modern version of the constituent quark model which is given 
by the coupled SD-BS approach, provides from = to — > oo the description for 
77* TT^ which is - independently of model details - consistent not only with the Abelian 
axial anomaly but also with the QCD predictions (^ for the leading large-Q^ behavior. In 
the SD-BS approach, is a quantity calculable through the straightforward application of 
the Mandelstam-formalism expression 

Pm/- = ^^/^^tr[x(g,p)7M75]- (11) 

The model dependence is present in the (successfully reproduced |]7|,|9|JT0|,|3[| ) value of but 
the derivation of the asymptotic forms (^ and (|T5|) below, is model- independent. They 
do not depend on what are the bound-state solutions, just like the chiral-limit axial-anomaly 
amplitude (|T]) doesn't. 



Of special importance is also that the derivation in Ref. [21| of the asymptotic large-Q 
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forms (^ and ([l5| ) below, applies to both Minkowski and Euclidean space. 



Our result ||2T| on the asymptotic behavior of Tt,()(—Q^, 0) subsequently received further 



support from Roberts |3^, who generalized its derivation by taking into account renormal- 



ization explicitly. His derivation shows that the asymptotics of Ref. [21] obtained using the 
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CP or mCP vertex, namely = 4/3, must be precisely reproduced for any gg7 vertex which 
is (like the CP or mCP ones) consistent with multiplicative renormalizability. 

It is interesting that the asymptotic behavior for large negative k"^ = —Q^ predicted ||21 



for the bare qq^f vertex and the dressed CP and mCP ones, is in exact agreement with the 
leading term predicted by Manohar using OPE [^]. In this paper, we shed some light on this 
connection by providing the following alternative derivation of the asymptotic behavior: con- 
sider TJ^(a;) J^(a;), the T-product of two quark vector currents (x) = ip{x)'y^{X°' /2)'ip{x), 
along the hues of Ref. ||3^ (Ch. 18 on OPE), but evaluated between the pion state and 
the vacuum. We substitute for the quark propagator its leading light cone [(x — y)^ — > 0] 
behavior 

S{x — y) ^ ^ \- — — , (12) 

[2n)^[[x — uY — lef 

and utilize 'ji^'f^'y'' = S^^^'^'y^ — is^^'^'^lal^ similarly as in the derivation of the large Q^- 
behavior in Ref. pl[. The term which contains the symmetric tensor S^^'^" = g^^g'^" + 



giJ-'^g^^'^ -)- g^^^g^" ^ is readily seen not to contribute to the 71*^77* vertex. On the other hand, 
in the term containing the antisymmetric Levi-Civita tensor £^^^'^" and 75, one finds - in 
the lowest order of the {x — ?/)-expansion - the pion-to- vacuum matrix element of the axial 
current, 

{m^)l,l^Q'mW\p)) = ^(0|^(0)7M75y^(0)|vr3(p)) = ^tUp,, (13) 
defining the pion decay constant The coordinate-space integration is reduced to 

d'-e^'-'7-I^.=^-'lJh:- (14) 



— ieY k"^ + 



Putting all this together reveals that the lowest-order result in the {x — ?/)-expansion is just 
Eq. @j with J being precisely 4/3. 

We thus in another way recover the result of Manohar pO| , of Ref. for the case when 



the bare (7^) or the dressed CP or mCP vertices are used, and of subsequent generalization 
| 33| thereof for any gg7 vertices consistent with the multiplicative renormalization. However, 



this asymptotic behavior Q'^Tt,o{—Q'^,Q) = (4/3) f-,, lies some 20% below the central values 
of the largest-Q^ CLEO data, and is more than one standard deviation away from them. 
Nevertheless, Manohar pointed out that his OPE approach also indicates the existence 



of potentially large corrections to his leading term. 

Let us now return to the BC vertex, which has been the most commonly used dressed qq^ 
vertex in the phenomenological SD-BS calculations so far. In conjunction with our SD-BS 
solutions, the usage of the BC vertex raises the asymptotic coefficient J' of Eq. (|[) from 4/3 
resulting from the usage of the bare, CP and mCP vertices, to ~ 1.78. The reason for 
this enhancement is that the BC vertex does not reduce to the bare vertex if the large 
momenta flow through just one, and not both, of its fermion legs. (Also note that the BC 
vertex is not consistent with multiplicative renormalizability p5|j29|] , so that the arguments 



of Ref. are not applicable to it.) It was shown in Ref. that when the BC vertex is 



used in our approach, the predicted asymptotic behavior is 
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T.o(-Q^O)-^A (15) 



where Z^- is given by the same expression (0) as /tt, except that the integrand is modified by 
the factor [1 + A{q^)Y' / A. In the special case of the solutions we use, the slow and moderate 
variation of A{q^), as well as its shape (see Fig. 3 in Ref. 0) permits the approximate 
factorization 

This last approximation (a rather rough one) would imply J ^ 1.69, i.e., 
Q2T^o(-g^0) ^ 0.157 GeV. The more accurate Eq. (0) gives Q'^T^o{-Q'^ 0.168 
GeV. This is practically indistinguishable from the asymptotics indicated by our full nu- 
merical calculation of the transition form factor, which reached Q'^T^^o^—Q'^ .,{]) 0.165 
GeV at the largest numerically studied . It is of special importance that analytically 
obtained large-Q^ asymptotics (|l^), the derivation of which is valid in both Minkowski and 
Euclidean space, makes this smooth and accurate contact between our large-Q^ numerical 
results calculated with the Euclidean solutions of the chosen model 0. The excellent agree- 
ment between the analytical and numerical results confirms the accuracy of our numerical 
methods and procedures (employing the Euclidean bound-state solutions 0) used in this 
and earlier papers i9|JlO|J^. 

Another case of this agreement between the numerical results and the analytical consid- 
erations is provided by (5^Tn-o(— Q^, 0) found in the soft and chiral limit approximation (^. 
Therefore, let us for a moment consider again this approximation, although we surpassed it 
already in Ref. |^ . (As a bonus, we also gain understanding of some subtle aspects in its pre- 



vious applications.) In Fig. 1, the results of the numerical calculation in that approximation, 
and with the BC vertex (|^), are depicted by the dash-dotted line, which practically reaches 
(from above) the asymptotic behavior, although not as fast as the results without that ap- 
proximation, depicted by the solid line. We now note that in the approximation (|^), the 
asymptotic behavior (|T3|) should be realized with and f-,, calculated in the same approxi- 
mation. We will denote them by f^^^ and f^^^, where gPS is short for "generalized Pagels- 
Stokar". Namely, it turns out that the soft limit approximation (|^), when applied to /^r. 



is the same as the widely used Pagels-Stokar approximation I^B] except that they work with 
the restriction A[q^) = 1 - unlike the present approach. Since we find f^^^ = 109.2 MeV, it 
appears that our dash-dotted curve should approach Q^T^f\-Q\0) (4/3)/^^^ = 0.146 
GeV. In fact, at = 35 GeV^, Q^Tf\-Q'^,0) ^ 0.136 GeV. To understand what forced 
the dash-dotted curve below, let us remember that the factor 1//^ appears in Eq. (0) be- 
cause in the chiral limit it turns out that the pion decay constant /^r is precisely equal to the 
the normalization of the BS vertex rT^{p,q) |3^. In Ref. the transition form factor was 
thus found using in Eq. (^ our chiral-limit value of the pion decay constant, /-n- = /° = 89.8 
MeV, and to make contact with that reference we have evaluated in this paper the dash- 
dotted curve in the same way. However, from what has just been explained, it follows that 
when we use not only the chiral, but also the soft limit (|^), the normalization of the BS 
vertex is given by the pion decay constant calculated in the same approximation, = f^^^, 
which in our adopted model has the value of 80.8 MeV. The transition form factor of Ref. P| 
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was thus suppressed by the factor f^^^/f^; for the same reason, the dash-dotted hne in the 
present Fig. 1 approaches from above (4/3)/^^'^(/^-^'^/ /°) ~ 0.131 GeV as the hmiting value, 
so everything talhes. 

In the case of the BC vertex, the "soft" leg adjacent to the pion BS-amplitude always 
contributes A{q'^), which enhances the asymptotic value of (5^T^o(— Q^, 0). This is not so for 
the CP or mCP vertices. Since they tend to the bare vertex as soon as the high momentum 
flows through one of their legs, does not get replaced by /tt, and the "bare" result, Eq. 
with J' = 4/3, continues to hold for Q"^ ^ oo when the CP or mCP vertices are used. 

Of course, the usage of the mCP vertex instead of the bare j'^ makes a considerable 
difference for the flnite Q^. Let us therefore now focus our attention to Fig. 3, which shows 
Q^Tt,o{—Q^, 0) found in our approach for both the BC-case (solid curve) and the mCP-case 
(dashed curve) for the interval of the transferred momenta between the presently highest 
accessible and = 0. 

At Q2 = 0, both BC and mCP vertices give the same amplitude (^ for 7r° — 77 in 
the chiral limit. With growing Q^, the BC- and mCP-curves in Fig. 3 must nevertheless 
ultimately separate because of the difference (4/3) (/,r — fn) of their large-Q^ limiting values. 
At the highest presently accessible momenta, ~ 8 GeV^, they differ by 9%, but as 
we go down in Q^, we see that already for < A GeV^, they practically coincide. This 
insensitivity on the choice of the WTI-preserving gg7-vertex is very indicative. Namely, the 
true solution for the r^{p,p') vertex, the one which we try to imitate by the BC and mCP 
vertex Ansdtze, is - for Q"^ < 4 GeV^ - unlikely to give (5^T^o(— Q^, 0) signiflcantly different 
from the BC and mCP vertices, since at = the transition amplitude must (in the chiral 
limit) again be ([l|), and at larger the difference should not be larger than that resulting 
from the BC and mCP vertices, as they lead to so very different Q^T^o(— Q^, 0) as — > 00. 
(One should also keep in mind that the multiplication by serves only for exposing the 
asymptotics more clearly, and that the difference is even smaller for Tj^o[—Q'^,0) proper, 
which must in every case decrease as l/Q"^.) The weak sensitivity (at least for not too large 
Q^) on the WTI-preserving Ansatz for r^{p,p') means that high-precision measurements of 
T^o(— (5^,0) will test unambiguously (and possibly give a hint on how to improve) the SD 
and BS model solutions which have so far been successful in fltting the low-energy hadron 
properties such as the meson spectrum. For example, our model choice |^ obviously (for 
both BC and mCP vertices) slightly overestimates T^o(— (^^,0) in the region < A GeV^. 

We point out that out of various SDK and BSE solutions in the SD-BS approach, once 
that /tt has been correctly reproduced by x [see Eq. (|11])], the transition form factor is most 
sensitive on v4(g^), or, more precisely, on its values at small and intermediate momenta — g^, 
where v4(g^) is still appreciably different from 1. Eqs. (P^-(P^ already explained how A[q'^) 
drives Q^T^o(— Q^, 0) upwards through /^r for large in the case of the BC vertex. 

To illustrate what happens at Q"^ < 8 GeV^ when the A(g^)-profile is decreased, let us 
enforce by hand the extreme, artiflcial case Alq"^) = 1. (To avoid confusion, we stress it is 
for illustrative purposes only, as we cannot have such a SD-solution in the adopted approach 
of Refs. |]23|j2^ ,[7[|.) This leads to the curve traced on Fig. 3 by small crosses, pertaining to 
the usage of both BC and mCP vertices (as well as the CP ones), since Alq"^) = 1 makes 
^mCP ~^ ^Bc- This curve illustrates how the heights of the curves depicting Q'^T^a{—Q'^, 0) 
depend on how much the 74(g^)-proflle exceeds 1. Obviously, for both the solid curve and the 
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dashed one, the agreement with experiment would be improved by lowering them somewhat 
(at least in the momentum region < 4 GeV), which could be achieved by modifying the 
model 1^ and/or its parameters so that such a new solution for A(g^) is somewhat lowered 
towards its asymptotic value oo) 1. (Of course, in order to be significant, this 

must not be a specialized re-fitting aimed only at A{q'^). Lowering of A{q'^) must be a result 
of a broad fit to many meson properties, comparable to the original fit 0. This, however, 
is beyond the scope of this paper.) 

By the same token, high precision measurements of T^o(— Q^, 0) (such as those planned 
at Jefferson Lab |]T3|]) can be especially helpful in obtaining information on A(g^) empirically. 
First note that Fig. 3 reveals that in the region around ~ 4 GeV^ we are already in the 
regime where the dependence on the BS-amplitude x is lumped into /,r- Therefore, high- 
precision measurements of T7ro(— Q^, 0) in the region close to the asymptotic regime will give 
information on the integrated strength of Alq"^) although not on A(g^) itself. However, since 
it is known that the form of that function must be a rather smooth transition {e.g., see [Q) 
from A(g^) > 1 for near 0, to ^(i?^) — 1 in the Q^-domain where QCD is perturbative, 
such measurements [|1^] would give a useful hint even about A{q'^) itself - namely about what 
solutions for v4(g^) one may have in sensible descriptions of dynamically dressed quarks and 
their bound states. 
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FIGURE CAPTIONS 

Fig. 1: Our results for (5^T^o(— Q^, 0) are presented from = up to = 36 GeV^. The 
data presently exist only for < 10 GeV^. CELLO data points are circles and those 
of CLEO are triangles. The dash-dotted line represents our Q'^T^o{—Q'^,0) evaluated 
in the chiral and soft limit approximation and only with the BC qq'j vertices. Our 
results without that approximation are depicted by the solid line for the case of the 
BC vertex, and by the dashed line for the case of the mCP vertex. This latter one gets 
lower beyond ~ 2.5 GeV^ and does not yet saturate, but at higher ultimately 
must approach asymptotically Q^T^o(— Q^, 0) — > 4/^/3 for higher Q^. This limit is also 
denoted by the dashed hne, but the straight one. The Brodsky-Lepage interpolation 
formula []T3| is the dotted curve approaching the (also dotted) straight line denoting 
the pQCD limiting value 2f.,^. The line of empty squares denotes basically the pQCD 
prediction (but with fixed Nf = 5) of Ref. [pT[]. 

Fig. 2: The pion distribution amplitudes (|TOD are depicted in the dimensionless version 0(a;) = 
fn{x) I for the three values of C pertinent to the discussion in this paper. C, = \ yields 
the asymptotic distribution amplitude (the dotted curve), C = 1-5 pertains to the pion 
distribution amplitude depicted by the solid curve and corresponds to the large-Q^ 
values of our T^o(— (^^,0) evaluated with the BC qq^ vertex. C, = 2.5 corresponds to 



the transition form factors predicted by Refs. [31] and [18| for the largest Q -values 



(Q^ ~ 8 GeV^) at which T^.o^—Q'^, 0) has been measured so far. 

Fig. 3: The comparison of our results for the pion transition form factor (times Q"^) with 
the CELLO (circles) and CLEO (triangles) data. Our results for Q^Tt^o[—Q'^,0) are 
depicted by the solid line for the case of the BC vertex, and by the dashed line for the 
case of the mCP vertex. This latter one gets lower beyond ~ 2.5 GeV^ and does not 
yet saturate at the presently accessible momenta although approaches asymptotically 
Q^T^o{—Q'^,0) — > 4/^/3 for much higher Q^. (This limit is denoted by the dashed 
straight line.) The little crosses denote (for both BC and mCP qq'j vertices) our 
(5^T^o(— Q^, 0) when we enforce A(g^) = 1 by hand. 
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